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Abstract. This thesis presents innovative contributions to the study of metal-dielectric- 
metal (MDM) integrated optical circuits with surface plasmon polaritons, using the 
transmission lines analogy for microwaves, the transfer matrix theory and the even-odd 
method. In the last chapter, I studied a plasmonic pulse to see how its shape changes 
when passing through waveguides with a variable number of periodic cells. 
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1. Introduction 


Surface plasmon polaritons (SPPs) are non-radiative electromagnetic waves, 
confined to the interface of a metal with a dielectric, obtained following the 
interaction of light with the collective oscillations of free electrons in the metal. 
Their applications are in many important fields such as data storage, Raman 
spectroscopy, materials with negative refractive index at frequencies in the visible 
region, microscopy, solar cell fabrication, waveguides and the detection of 
biomolecules using sensors [1]. 


Integrated optical circuits could use SPPs for all the circuit’s elements, which 
results in a miniaturization of circuits under the diffraction limit. SPPs can be 
used, for example, to fabricate very compact nanophotonic circuits and for optical 
interconnection on a motherboard. The study of integrated optical circuits with 
SPPs has many applications, which offer the possibility of theoretical 
investigation and practical realization of optical devices that can operate at the 
nanometric level [2]. 


2. Surface Plasmon Polaritons 


Monochromatic transverse magnetic plasmonic waves (with frequency @) are 
obtained at a metal-dielectric interface, as shown in Figure 1. This figure also 
illustrates the spatial distribution of the electric and magnetic fields of the SPP. 


In agreement with the condition of continuity at the interface for the tangential 
components of the electric field, we find the dispersion relation of SPPs: 
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Fig. 1. Representation of the distribution of the components of the electromagnetic field of a 
plasmonic wave that appears at an interface between a metal and a dielectric 
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If we graphically represent the dispersion relation (1) (see Figure 2), three areas of 
interest are distinguished: in the first domain @ < w@spp for which ém < -éa 1.e. for 
frequencies w lower than the resonance frequency of the surface plasmons qwspp, 
we have kspp > ko and the dispersion curve - represented by a black line, is below 
the curve corresponding to the incident radiation, i.e. to the curve ko = a/c - 
marked with a black dotted line on the graph. The electromagnetic field of the 
SPPs decreases exponentially in the normal direction at the interface. The second 
domain is the one for which wspp < w < wp, where wp is the plasma frequency. 
Here, in the absence of any loss, kspp has only an imaginary component, which 
means that the plasmonic wave does not propagate along the interface. The third 
domain is the one for which @ > wp, where we have volume plasmons. In the 
simulation we took é¢ = éair = 1, and 29 = 0.632 wm the wavelength of the 
electromagnetic radiation in vacuum. Regarding the permittivity of silver, this is 
described by the Drude-Lorentz model [3]. All simulations in this thesis were 
done with the Scilab open-source program. In all the simulations performed in this 
thesis, the losses were neglected and only the real part of the electrical 
permittivity dependence of frequency was considered: 


2 (a) =1-—, (1) 


where wp = ne’/eom is the plasma frequency. In the case of silver, &m = Esiver = 1- 
wp'/a@* with wp = 8.85x10'> Hz [3]. 
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Fig. 2. Representation of the dispersion relation of the SPP for an air-silver interface 


The propagation length of the SPP is given by: 
Lspp =1/Im(2k spp ) (3) 


In practice, losses in conductors, called ohmic losses, limit the propagation length, 
which is a problem for the components of the SPP waveguides. The attenuation 
distances in dielectric and metal (for an attenuation of 1/e from the initial 
amplitude) are given by the relations: 


A | €,+6é A E, +é 
d = m d d —_—~ {_“m d 4 
d on e m nx E ( ) 


=>e,-d, =|€,,|-d,,- Because | ¢,,] >>, , we can say that d,, <<d,. 


m 


In practice, dm is about tens of nm and weakly depends on the wavelength J, 
whereas da grows rapidly and nonlinearly with /: 


d, ~cla, d, =(4/2)’ -(0, Ice,) (5) 
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SPPs cannot be produced directly because, for the same frequency, a phase 
mismatch occurs between the incident light vector ko and that of the plasmon- 
polariton denoted by kspp. Because we have kspp > ko, this translates to the fact 
that SPPs cannot be excited by direct metal illumination. Various methods have 
been developed to overcome this mismatch, like Kretschmann configuration with 
a prism, Kretschmann configuration with a two-layer prism, difraction grating, 
microscope probe with scanning in near field SNOM or Otto configuration with a 
prism [1]. 


3. The Analogy with the Microwave Transmission Lines 


MDM (metal-dielectric-metal) and DMD (dielectric-metal-dielectric) are the most 
common types of plasmonic waveguides, each having its characteristics, 
advantages and disadvantages [5]. Thus, in a MDM type waveguide, the magnetic 
field is more intense in the dielectric region and weaker in the metal region, 
because the plasmonic waves that appear at the two metal-dielectric interfaces 
overlap when the width of the guide is small, while in a DMD waveguide, we 
have the opposite situation. Thus, the magnetic field is located almost entirely 
inside the guide, with a minimum penetration depth outside the cross section of 
the MDM waveguide. In addition, the propagation length in MDM can exceed 
tens of microns for waveguides with widths less than 10% of the wavelength of 
vacuum light [4]. For DMD structures, the propagation length of the SPPs is of 
the order of centimeters for wavelengths in the near IR range, but the penetration 
depth of the SPPs is of the order of microns in the dielectric zone. [5] 


Frequently MDM plasmonic waveguides are described by solving Maxwell 
equations by FDTD (finite difference time — domain) method, but this is time 
consuming. Another way to model SPP propagation through metal-dielectric- 
metal waveguides is to use the analogy with the microwave transmission lines. In 
this respect, it is possible to propose equivalent circuits for the MDM waveguides 
and to calculate the transmission through the waveguide, using the equations of 
the transmission lines. 


In this thesis I studied various MDM configurations, including a single-arm MDM 
plasmonic waveguide, two MDM guides connected by a coupling section and an 
architecture consisting of two MDM _ waveguides, with an arm coupled 
perpendicular to the first guide. For example, the case of a single-arm MDM 
plasmonic waveguide is represented in figure 3, which indicates also the 
geometric parameters of the configuration. 


If the dimensions of the waveguide and the arm comply with the conditions h << 4 
and respectively w << A, only one transversal TM mode will propagate through 
the waveguide. 
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Fig. 3. (a) Schematic of a MDM waveguide with width h, with a single irregularity coupled 
perpendicular to the guide, length d and width w. (b) Representation of the circuit equivalent to (a) 
according to the transmission line theory and (c) its simplified representation 


We can schematically represent the waveguide studied with a network consisting 
of an infinite transmission line with characteristic impedance Zo coupled in series 
with a finite transmission line with characteristic impedance Zarm and we can write 
the transmittance at the end of the MDM guide using the transfer matrix analysis 
method, as the inverse of the square module of the element (1,1) of the total 
transfer matrix, T= 1/Mii/’: 


Transmittance 
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Fig. 4. Transmission curves for two Ag-air-Ag waveguides, one with an arm coupled 
perpendicular to the guide with ¢ = 1 (continuous line) as in the rest of the guide, and another with 
the arm with ¢ = 1.2 different from the rest of the guide (with dotted line) 
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where Ai(/) is the matrix representing the propagation of SPPs through the 
waveguide with length /, B is the matrix representing the reflection of SPPs on the 
non-uniform section coupled perpendicular to the first guide, and A:(L-/) is the 
matrix that models the propagation of SPPs through the guide of length L-I. 


I compared two transmission curves, one for an Ag-air-Ag waveguide with one 
arm coupled perpendicular to the waveguide with ¢ = | as in the rest of the guide, 
and another for a waveguide with one arm coupled perpendicular to guide with ¢ = 
1.2 (see figure 4). From simulations, it turned out that the electrical permittivity 
and the arm length are the parameters that influence the position of the maximum 
transmittance, that is, we have a displacement of the curve towards smaller 
wavelengths when the electrical permittivity in the arm increases compared to the 
rest of the guide or, respectively, the arm length increases. 


4. Plasmonic Logic Gates 


For the implementation of some logic gates a MDM slot waveguide with constant 
section and one, two or three regions with variable electrical permittivity or 
different refractive index can be used (see Figure 5, which also indicates the 
dimensions of the configuration) [8]. Graphene can be embedded in the slot 
waveguide to obtain a variable permittivity region. This structure is 
reconfigurable. Reconfiguration involves applying another set of gate voltages on 
the three regions, without the need to change the dimensions of the guide. 
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Fig. 5. MDM waveguide with three sections covered with graphene 


For example, I simulated a three-bit Toffoli reversible universal logic gate using 
this structure. In this case, the first graphene region represents the first bit, with 
the logical input states 0 and 1, given by refractive indices corresponding to &go 
and &¢1;, while the following two graphene sections represent control bits with 
values 0 and 1, encoded by the values of the refractive indices corresponding to 
Eedc and gic. I represented in the figure 6 the transmission coefficient as a function 
of the wavelength, for the following parameters of the structure: /2 = /3 = /g1 = 
lg2 = 1g3 = 100 nm, d1 = d2 = d3 = 40 nm, dg = 40 nm, ég = 1.2, &¢1:= 5 when 
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the first bit is O and égoc = 4.75, égic = 1.2, that is, we represented states 000 (with 
black dashed line), 010 (with dash-dotted line), 001 (with dash-circle line) and 
011 (with solid black line). 


Transmittance 
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Fig. 6. Simulation of the Toffoli reversible universal logic gate for states 011 (with black dashed 
line), 010 (with dash-dotted line), 001 (with dash-circle line) and 000 (with solid black line) 
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Fig. 7. Simulation of the Toffoli reversible universal logic gate for states 100 (+ black line), 110 
(black dash-circle line), 101 (> black line) and 111 (black solid line) 


In Figure 7, I represented the transmission coefficient as a function of the 
wavelength, for the following structure parameters: /1 = /2 = 13 = /4=/g1 =1g2 = 
lg3 = 100 nm, d1 = d2 = d3 = 40 nm, dg = 40 nm, é¢0; = 1.2, &g1+= 5 when the first 
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bit is 1 and €goc = 1.2, Egic = 4.75, that is, I represented states 100 (+ black line), 
110 (black dash-circle line), 101 (> black line) and 111 (black solid line). 


From these two figures it can be understand the operating principle of the Toffoli 
gate, which inverts the value of a bit when both control bits are 0. For this to 
happen, the value of transmittance T must have different values for 0 and | states. 
In this situation T is 0 for low transmittance, i.e. for T < 0.7 and 1 for T > 0.9 in 
Figure 6 and has opposite meaning in Figure 7. The transmittance changes from 0 
to 1 state when the last bits are 11 in Figure 6, i.e. Tis small for 000, 010 and 001 
states and large for 011 state, and from | to 0 in Figure 7 if the last two bits are 
11. The logic gate in these figures operates at 1.55 um, which is an important 
wavelength in applications. Other logic gates can be implemented, for example 
CNOT, AND, OR if two graphene sections are activated by gate electrodes [8]. 


5. Ring Shaped Plasmonic Logic Circuits 


The three- and four-port ring structures could function as various reversible and 
irreversible logic gates, if the amplitude and phase of the input signals are 
changed. The ring-type configurations studied in this thesis are symmetric to be 
able to use the even-odd analysis method of transmission lines [9]. For example, I 
investigated logical gates implemented through a symmetrical configuration with 
three ports like the one in Figure 8 [12]. 


Port 3 


Port | 


Port 2 


Fig. 8. MDM ring-shaped plasmonic waveguide, with three ports 


a. Schematic representation b. Even-odd analysis equivalent circuit 


The ring was chosen so that /2 = 2/; and the input signals at ports 1 and 2 are a; 
and azexp¢, where aj, a2 are the amplitudes of the signals, ¢ is the relative phase 
between the two inputs, and the complex output signal from port 3 is denoted as t. 
The transmittance is T = |t\’. The slot waveguide has a structure Ag/air/Ag, the 
relative dielectric constant of air is considered equal to 1, while the dielectric 
constant of the silver ém, is taken from [3]. In the following simulations was 
considered a waveguide with a width d = 50 nm, illuminated at Aexe = 550 nm, 
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with kspp = 1.68x10’ m'! greater than kexe = 27/ deve = 1.14x107 m'!; the matching 
of the wave numbers and, implicitly, the excitation of SPPs can be done, for 
example, by means of a diffraction grating. In this case, we can replace the slot 
waveguide with a transmission line [6-7] with characteristic impedance given by: 


2 llegeaid 


Z.= 
K eCEq Eo 


= 27,6 (7) 


0 


The transmittance T is a periodic function of the length /; (see Figure 9) for ai = 1, 
a2 = 1, ¢ = 0 (solid line), a; = 1, a2 = 0 or a; = 0, a2 = 1, regardless of ¢ (dashed 
line) and a) = 1, a2 = 1, 6=7 or a; = 0, az = 0 (dashed-dotted line). 
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Fig. 9. Dependence of T on the dimensions of a 3-port ring circuit 


Each time [jkspp is an integer multiple of 2, the transmittance is maximum. This 
maximum is four times higher than in the case of a single input signal if the two 
inputs are in phase or disappear when the two inputs are out of phase (see Figure 
9). This result can be analyzed by analogy with a ring-type acoustic gate [10], but 
in this case we measure intensities, and not acoustic amplitudes so that the 
variation of the signal from port 3 is much larger than in the acoustic case. 


For a maximum 7, suppose that /ikspp is an integer multiple of 2, and logical states 
0 and 1 are encoded in the amplitudes of input states at ports 1 and 2. Then, as for 
the acoustic logic gate [10], the 3-port circuit implements the irreversible gate OR 
for ¢ = 0 (the inputs in phase) if the output state is associated to logic values as 
follows: logic 1 for T > 0.3 and O otherwise. Indeed, in this case, as shown in 
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Figure 9, at output we have logic state | if one or both inputs are | logic and 0 if 
both are 0 logic. But there is a potential problem with the reading of output state: 
the transmittance has the same logic value for T about 1.8 when (a1, a2) = (1,1), 
and about 0.45 when (a1, a2) = (0,1) or (1,0). This very high difference between 
these output transmittances associated to logic state 1 can be attenuated through a 
different ¢ value. The phase difference ¢ can be changed by modifying the lengths 
of the input ports with AL = L2 - Li = d/kspp. 


Similar, for ¢ = a (the inputs at ports | and 2 are out-of-phase), Figure 9 suggests 
that the circuit implements the XOR gate, which is an irreversible logic gate, for 
the same designation of the output state at port 3: T is associated to logic state | 
only for inputs (a1, a2) = (0,1) or (1,0). 


If in figure 9 we associate the logic state 1 with T > 1 and the logic state 0 with T 
< 1, the circuit works as the gate AND, which is again an irreversible gate, for ¢ = 
0. In this situation the output state is 1 only when (a1, a2) = (1,1). Also, the 
configuration of the three-port slot waveguide studied above can implement the 
CNOT reversible logic gate. 


I also studied the output intensity at port 4 of a four-port ring-shaped symmetrical 
waveguide, finding that the dependence of the output signal on the circuit's 
dimensions is more sensitive to their changes for the 4-port configuration than for 
a 3-port configuration. This is a situation that is not desirable because it imposes 
strict manufacturing tolerances. To avoid it, the 4-port circuit could be replaced by 
two 3-port ring-shaped plasmonic waveguides placed in cascade, for which the 
tolerances are more flexible [12]. 


6. Plasmonic Pulse Shaping 


In this chapter is shown that the shape of a plasmonic pulse can also be changed 
using plasmonic circuits which consist of slot waveguides and arms, not only, as 
usual, by the controlled modulation of the spectral components of the optical 
pulse, which are spatially separated by gratings [11]. 


I studied a MDM slot waveguide, with the dielectric region having the width d 
and the electrical permittivity ¢a, and the metallic regions with the electrical 
permittivity ém. The waveguide has different arms attached in a periodic cell 
structure, the configuration of a single cell being represented in figure 10. The 
plasmonic circuit has N cells, each with two arms with the lengths /; and Jz 
respectively at the distance L, and the widths are denoted by d; and d2 
respectively. The distance between the center of the second arm and the end of the 
cell is /, so that the distance between two groups of arms in adjacent cells is AL =/ 
+ di/2. 
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Fig. 10. The structure of a single cell 


In order to obtain a much different shaped output pulse from that of the input 
pulse, the optical circuit must have a strongly frequency dependent transmission 
coefficient and the incident optical pulse should be as short as possible. I used the 
analogy with the microwave transmission lines [6-7] to describe each section of 
the waveguide length L through a matrix A(L): 


exp(—ik ppl 0 
A(L) = P(-tk sppL) . (8) 
0 exp(ik sppL) 
and the impedance given by 
For each arm /, where j = 1, 2, ... the matrix Mj can be written as follows: 
1+Z,,;/2Z, =Z,;/2Z, 
A Oh, Waa A OZ 
Z, =K spp jd ,/@&, aes 
®;=Kspp, jf 


VE, /E, —itang, 
Z,5 =4) 
l-ijé,/é,, tang, 
where kspp, ; is the plasmon propagation constant through the arm width dj, with j 


= 1.2. Also, I considered that t(a) =1/M,,. where M1: is the element (1, 1) of the 
matrix that describes the plasmonic circuit, namely the matrix 


M=[M , A(L)M , A(AL)]* (11) 
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The simulations showed that the shape of the incident pulse can be changed by 
changing the dimensions of the circuit for a slot waveguide, with an Ag-air-Ag 
configuration, through which is transmitted a Gaussian input pulse f(t) = exp(- 
?/2Ar’) with pulse duration At and the central frequency wo, corresponding to the 
wavelength of 800 nm and the width of the waveguide d = 70 nm. The circuit 
consisting of N cells modifies the pulse based on the following principle: The 
Fourier transform of the input pulse is multiplies by t(@), so that the output pulse, 
obtained by taking an inverse Fourier transform, differs from the incident pulse. 


For example, the left column of Figure 11 indicates with black dashed line the 
absolute value of the incident pulse, normalized to its maximum value, and with 
continuous black line the absolute value of the output pulse, while in the right 
column there are, the absolute values of the Fourier transforms for the input pulse 
(with black dashed lines) and the output pulse (with continuous black lines), 
normalized to the maximum of the Fourier transform of the incident pulse. The 
transmittance T = |¢|? (blue line) is also indicated. This figure shows that the output 
pulse is shorter than the input pulse and has an asymmetric shape because the 
transmittance has a strong variation with frequency around wo. The output pulse is 
smaller than the incident pulse because of light reflection. 
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Fig. 11. Pulse evolution through the plasmonic circuit with At = 50s, N = 3, d;) = 75 nm, d2= 50 

nm, /; = 450 nm, 2 = 315 nm, L= 150 nm and AL = 200 nm 


To study the influence of each parameter of the plasmonic circuit on the 
transmittance and, automatically, on the output pulse shape, I varied dj, d2, hi, lo, 
L, AL in turn and observed the effect of these variations on T and the shape of the 
pulse. From simulations it was observed that the most important parameter to 
modify if we want to change the output pulse is the length of the arms. For 
example, if we want to obtain a displacement with Aw of the transmittance, we 
can, alternatively, vary by about 7% the distance between the arms, by about 9% 
the width of the arms or by about 0.4% their length. The applications of pulse 
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shaping at the exit of a plasmonic circuit range from microscopy, imaging, 
spectroscopy, high-speed communications to compressing optical pulses [18-23]. 


Conclusions 


Surface plasmon-integrated optical circuits represent an alternative to electronic 
circuits [13-17] and have applications in fields such as waveguides, microscopy, 
data storage, Raman spectroscopy and biomolecules detection. Even though they 
are not as widespread today as electronic circuits, optical and especially 
plasmonic circuits are of interest because they can transmit information at higher 
speeds and, they are very compact. 


The novelty of this thesis is the demonstration that architectures based on MDM 
slot waveguides can be used to implement logic gates, if they are excited by 
monochromatic radiation, or they can change the shape of light pulses. The 
simulations that prove these functionalities were based on the analogy with the 
microwave transmission lines and were performed using the open-source software 
Scilab. 


It was shown that logic gates of the type NOT, OR, AND, CNOT and Toffoli can 
be implemented in a MDM plasmonic waveguide, with one, two or three sections 
covered by graphene. All of these logic gates are reversible and can be 
implemented in a single configuration containing three graphene regions whose 
refractive index can be adjusted independently using individually contacted gate 
electrodes. The results obtained in this chapter have been published in [8]. In all 
situations we have found that transmission can be strongly modulated, which is 
important for practical applications involving logic gates. 


Innovative is also the use of analogy with the microwave transmision lines to 
demonstrate that ring-shaped, three- and four-port plasmonic slot waveguides can 
implement multiple logic gates with one, two and three bits depending on their 
configuration, amplitudes and phases of the input signals. These results have been 
published in [12]. In particular, we have shown that two ring structures with three 
ports can be cascaded in order to implement the universal and reversible CCNOT 
logic gate in a more efficient way than a similar ring structure, but with four ports. 


Another original approach of this thesis concerns the demonstration of the 
possibility of adjusting the shape and the duration of an incident optical pulse 
when passing through plasmonic circuits consisting of periodic structures of arms. 
These changes are based on the different modulation of the amplitude and phase 
of the spectral components of the input pulse. We showed that the shape of the 
output pulse is influenced most by the variation of the arm length. These results 
are sent for publication. The particular importance of obtaining optical pulses with 
different shapes consists of possible applications in fields such as microscopy, 
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optical pulse compression, imaging, high speed communications and 
spectroscopy. 
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